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摘 要 : 干旱 作为 频 发 的 全 球 性 自然 灾害 之 一 ,造成 了 严重 的 社会 、 经 济 和 生态 环境 问题 。 以 西藏 


主要 耕作 


区 为 研究 


区 ,2001 一 2015 年 MODIS、TRMM 和 SRTM-DEM 数据 为 数据 源 ,利用 植被 状态 


指数 (VCI) 、 温 度 条 件 指 数 (TCI1) 和 降水 状态 指数 (PC1) 等 模型 参量 ,采用 空间 主 成 分 分 析 方 法 构 


建 


部 及 东部 地 区 ; 


区 干旱 易 发 月 份 存在 差异 ; 
现 出 随 月 份 变化 由 耕作 
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干旱 作为 频 发 的 全 球 性 自然 灾害 之 一 ,造成 了 
严重 的 社会 .经 济 和 生态 环境 问题 ,对 区 域 粮食 安全 
和 社会 稳定 构成 严重 威胁 '"。 我 国 作为 农业 大 国 ， 
长 年 以 来 饱 受 干旱 影响 。 青 藏 高 原 地 区 农业 生产 及 
作物 生长 与 区 域 复杂 的 气候 环境 、 地 形 构 造 紧 密 相 
关 。 因 此 ,针对 高 原 地 区 独特 的 地 理 位 置 与 气候 环 
境 ,研究 适合 的 旱情 监测 模型 显得 尤为 必要 。 

由 于 干旱 受到 多 种 因素 影响 ,单一 要 素 往往 难 
以 表征 全 面 的 农业 干旱 信息 。 近 年 来 ,很 多 学 者 尝 
试 利用 多 源 遥 感 数据 构建 全 面 的 干旱 指数 模型 来 进 
一 步 提高 干旱 监测 精度 。KOGAN 4§°) 基于 多 年 的 
植被 指数 数据 ,通过 提出 植被 状态 指数 (VC7) ,并 结 
合 VCI 与 温度 条 件 指数 (7TC1) 提出 更 具 优势 的 植被 
健康 指数 (VHI) 来 进行 旱情 监测 。TADESSE 等 中 
利用 分 类 回归 树 模型 (CART) 数据 挖掘 的 方法 结合 
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区 域 干 旱 综 合 监 测 模型 ,对 模型 精度 和 可 靠 性 进 
2015 年 逐 月 干旱 进行 识别 ,采用 地 理 时 空 分 析 方 法 对 研究 
表明 ,干旱 综合 监测 指数 (DCMI) 能 够 较 好 地 反映 
(SPED) 的 变化 ,干旱 综合 监测 模型 具有 较 好 的 适用 性 ;研究 
西高 东 低 的 特征 ,大 部 分 地 区 干旱 频率 小 于 20% , 约 12.4195 的 
级 干旱 发 生 频 率 来 看 ,日 喀 则 市 为 轻 旱 .中 旱 易 发 区 , 重 旱 易 发 区 则 集中 于 日 喀 则 市 和 昌都 市 的 中 
区 域 月 际 干旱 频率 空间 格局 差异 较 大 ,全 年 干旱 易 发 生 于 1、8、11 月 等 月 份 , 局 部 地 
区 域 年 内 旱情 变化 趋势 差异 性 较 大 ,10 月 ~ 次 年 9 月 ,旱情 加 剧 区 
区 东部 向 西部 逐渐 转移 的 趋势 。 


1000 -6060(2019 )04 -0735 -10(0735 ~ 0744) 


中 国 科学 院 战 略 性 先导 科技 专项 (XDA20030302 ) ;水 科 院 全 国 山洪 灾害 调查 评价 项 目 
藏 自治 区 科技 支撑 计划 项 目 (省 809) ; 

了 信息 系统 理论 与 灾害 风险 分 析 方 面 的 研究 . E-mail; neu_xjn@ 163. com 
程 维 明 (1973 - ) , 男 ,人 研究员, 主要 从 事 数 字 地 形 地 貌 研 究 . E-mail; chengwm@ lreis. ac. cn 


行 检 核验 证 ,并 以 所 建 模型 对 研究 区 2001 一 
区 干旱 交 化 特征 及 趋势 进行 研究 。 结 果 
区 域 土壤 相对 湿度 与 标准 化 降水 茹 散 指数 
区 全 年 干旱 频率 在 空间 分 布 上 呈现 出 
区 域 干 旱 频 率 超 过 20% ;从 不 同等 


遥感 数据 来 进行 干旱 监测 。RAJSEKHAR 等 中 综合 
考虑 降水 .蒸发 以 及 土壤 湿度 等 因素 提出 多 变量 干 
旱 指数 (MD7) 用 于 气象 干旱 .农业 干旱 以 及 水 文王 
旱 的 监测 。ZHANG 等 钻研 究 发 现 与 其 他 单一 指数 
相 比 ,降水 状态 指数 (PC7) 几乎 与 所 有 短期 干旱 的 
相关 性 更 强 。 不 少 研究 还 发 现 ,土地 覆盖 类 型 (ZC) 
是 土壤 水 分 空间 变异 的 主要 影响 因子 '9; 使 用 
TRMM MODIS 等 数据 结合 气象 站 点 干旱 指数 ,利用 
数据 挖掘 算法 可 以 较 好 地 建立 综合 干旱 监测 模 
型 了 -8 。 此 外 , 杜 灵 通 等 "结合 降水 .地表 温度 和 植 
被 指数 数据 构建 出 综合 干旱 监测 指数 (SD7) ,并 取 
得 较 好 的 监测 效果 。 然 而 ,目前 国内 针对 区 域 干旱 
监测 的 研究 大 多 集中 于 干旱 监测 指数 、 监 测 模型 构 
建 与 验证 ,干旱 发 生 过 程 及 影响 范围 的 探讨 ,而 采用 
多 源 长 时 序数 据 建立 综合 干旱 监测 模型 ,对 研究 区 
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长 时 间 周 期 内 逐 月 干旱 进行 识别 、 分 析 其 时 空 变化 
规律 及 趋势 等 方面 的 研究 还 相对 较 少 。 而 历史 干旱 
识别 .时 空格 局 及 变化 趋势 分 析 ,是 区 域 干 旱 研究 与 
应 对 的 重要 内 容 。 

选取 地 处 西藏 自治 区 东南 部 的 拉萨 日喀则 
等 5 市 为 研究 区 ,以 MODIS、TRMM 和 SRTM-DEM 
为 数据 源 , 综 合 分 析 人 研究 区 气象 降水 亏 缺 .地 表 蒜 
散 、 地 形 等 因素 ,选取 VCI TCI 和 PCI 等 模型 参 
量 ,采用 空间 主 成 分 分 析 方 法 构建 了 人 研究 区 干旱 
综合 监测 模型 ,并 对 模型 的 精度 和 可 靠 性 进行 检 
核验 证 。 在 此 基础 上 对 研究 区 2001 一 2015 年 逐 
月 干旱 进行 识别 和 重建 。 通 过 时 空 分析 方 法 ,分 
析 研 究 区 历史 干旱 时 空格 局 及 其 演化 趋势 , 旨 在 
为 青藏 高 原 地 区 历史 干旱 时 空 规律 与 区 域 干旱 监 
测 预测 的 研究 提供 借鉴 ,为 开展 区 域 防 旱 抗旱 工 
作 提 供 科学 参考 。 


1 研究 区 概况 及 数据 源 


1.1 研究 区 概况 

研究 区 位 于 青藏 高 原 的 东南 部 .西藏 自治 区 
的 藏 南 谷 地 和 藏 东 峡 谷地 ,主要 包括 拉萨 .昌都 、 
林芝 山南 和 日 喀 则 5 个 市 (图 1) 的 56 个 区 
(县 ) ,区 域 面积 约 51.69 x 10° km ,平均 海拔 均 在 
3 000 m 以 上 ,耕种 土壤 面积 占 全 区 耕种 土壤 面积 
的 97.43% ,是 西藏 自治 区 最 主要 的 农业 区 。 人 研究 
区 内 气候 分 布 多 样 ,光照 充足 ,年 均 气 温 3.6 C, 
年 均 降 水 量 420 mm 左右 ,全 年 日 照 时 间 均 在 
2 000 h 以 上 ,无 霜 期 100 ~ 180 d。 据 不 完全 统计 ， 
1983—2013 年 间 全 区 共 发 生 大 小 旱灾 350 RK, 
其 中 拉萨 市 .昌都 市 等 5 市 共计 308 w, BEKE 
灾 总 数 的 85% 以 上 。 
1.2 数据 源 与 预 处 理 

研究 数据 包括 覆盖 研究 区 的 2001—2015 年 的 
植被 指数 数据 (MOD13A3 )、 地 表 温 度数 据 
(MOD11A2) .土地 覆盖 数据 (MCD12Q1 ) 、 降 水 数据 
(TRMM3B43) .SRTM-DEM 数据 和 气象 站 点 数据 ， 
数据 主要 来 源 如 表 1 所 示 。 由 于 MCD12Q1 HW% 
盖 数 据 只 更 新 至 2013 年 , 故 文 中 2014 年 .2015 年 
WRA tt BE A 2013 年 替代 。 

人 研究 过 程 中 利用 MRT ( MODIS Reprojection 
Tool ) 软件 对 MODIS 数据 进行 拼接 .投影 转换 ,应 用 
ArcGIS10.2 对 影像 数据 进行 裁剪 ,为 保证 数据 的 
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Fig.1 Major rivers and elevations in the study area 


表 1 监测 模型 主要 数据 来 源 
Tab.1 Monitoring model main data source 


空间 ”时间 


È 1 
数据 名 称 分辩 a A 数据 来 源 

植被 指数 产品 ik 月 值 5 景 / 期 x 12 https://mirador. 

(MOD13A3 ) /年 x15 a gsfc. nasa. gov 
地 表 温 度 产品 lk gd 5 景 /期 x 46 https://mirador. 

(MOD11A2 ) m /年 x15 a gsfc. nasa. gov 
土地 覆盖 产品 _ https://lpdaac. 

A 一 : 3 年 
(MCD12Q1) 500m 年 值 2001—201 ese 


降水 产品 https ;// mirador. 
5° 12 期 / 15 : 
(TRMM3B43 ) 0.25” 月 什 AVE x15 a gsfc. nasa. gov 
数字 高 程 模型 ik 1 http://www. 
(SRTM-DEM) ii cgiar-csi. org 
气象 站 点 2 £ 1 http :// www. resdc. 


cn/ 


完整 性 ,使 用 质量 控制 文件 剔除 MODIS 影像 中 的 无 
效 值 , 并 利用 多 年 均值 填充 算法 对 同 区 域 无 效 值 进 
行 补充 。 为 统一 数据 时 空 尺 度 ,将 8 d 地 表 温 度数 
据 ,利用 最 大 值 合 成 法 生成 地 表 温 度 月 值 数据 ,将 地 
表 履 盖 产 品 重 采样 为 1 km 分 辨 率 。 利 用 ENVI5.1 
对 TRMM3B43 进行 预 处 理 ,并 结合 增强 型 植被 指 
Bc" (EVI) RH CWR 模型 对 合成 月 降水 数据 进行 
降 尺 度 处 理 。 


2 干旱 监测 模型 构建 及 验证 


综合 考虑 研究 区 植被 覆盖 状况 .地 表 温 度 .降水 
情况 .土地 覆盖 以 及 高 程 分 布 5 种 影响 因子 ,本文 将 
X VCI TCI, PCI, LC, DEM 来 构建 高 原 地 区 干旱 
指数 监测 模型 ,模型 参量 如 表 2 所 示 。 
2.1 模型 参量 计算 
2.1.1 VCI 指数 植被 通常 对 干旱 有 着 明显 的 响 
应 。 当 植被 叶片 枯萎 ,凋零 和 和 履 盖 减少 时 均 会 引起 
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表 2 ”监测 模型 构建 指标 选取 
Tab.2 Monitoring model construction index selection 
干旱 影 指标 
响 因 素 选取 
植被 VCI ”该 指数 是 衡量 植被 受 环境 胁迫 程度 的 指标 , 反 
覆盖 映 短期 天 气 信号 对 植被 的 影响 ,能 有 效 监测 干 
旱 及 降水 的 时 空 分 布 特征 
气候。 TCI ”该 指数 强调 温度 条 件 与 植被 长 势 的 关系 ,在 植 
状况 被 稀 玻 及 作物 播 收 期 均 可 监测 , 视 为 旱情 监测 
的 初始 指标 


指标 含义 


PCL ”干旱 发 生 的 直接 原因 大 部 分 为 降水 异常 ,该 指 
数 可 有 效 监测 干旱 及 降水 的 时 空 分 布 特征 
地 表 LC KE LC 类 型 会 改变 土壤 性 质 及 土壤 结构 ,影响 
环境 土壤 的 持 水 能 力 和 水 分 调节 能 力 ,从 而 影响 十 
SES TCL 
地 势 DEM ”特殊 的 地 形 可 形成 独特 局 部 的 小 气候 ,间接 影 
影响 响 土壤 水 分 的 含量 和 分 布 
植被 敏感 波段 反射 率 的 变化 ,从 而 可 以 利用 听 感 影 


像 进 行 探 测 。 增 强 型 植被 指数 (EVI) 与 归 一 化 植被 
指数 (NDVI) 相 比 ,具有 比 NDVI 更 高 的 稳定 性 ,能 
更 有 效 的 反映 植被 覆盖 与 变化 情况 。 由 于 EVI 
难以 反映 时 间 序 列 上 植被 的 长 势 变 化 ,因此 本 文选 
择 由 KOGAN!) 提出 的 VCI 来 反映 植被 的 生长 状 
况 , 其 计算 公式 为 "|. 


VCI EVI, E EVI nin 1 
ar EVI pax = EVI rin | ) 


IP: VCL 为 上 月 的 植被 状态 指数 ;ETV7 为 i 月 的 增 
强 型 植被 指数 值 ;EV 为 多 年 时 间 序 列 中 i 月 增强 
型 植被 指数 的 最 大 值 ;EVL,,,, 为 多 年 时 间 序 列 中 i 月 
增强 型 植被 指数 的 最 小 值 。 其 值 域 范围 在 [0,1] 之 
间 ,通常 高 值 表示 该 月 较为 湿润 ,植被 生长 良好 , 反 
之 亦 然 。 

2.1.2 TCI 指数 已 有 人 研究 表明 ,干旱 通常 会 伴随 
着 明显 的 植被 冠 层 温度 异常 ”|。7CI 可 通过 遥感 
的 热 应 力 指标 地 表 温 度 (ZS7) 来 确定 相关 的 干旱 现 
象 。 其 原理 是 植被 冠 层 或 土壤 表面 温度 随 着 水 分 胁 
迫 的 增加 而 增加 ,不 受 作物 生长 季节 的 限制 ,适用 于 
长 时 间 序 列 及 大 区 域 的 相对 干旱 监测 "1 BD 
LST, — LST, 

LST — LST nin 


式 中 :7C1 为 了 月 的 温度 条 件 指数 ;2S7 为 i 月 的 地 
表 温 度 值 ;LS7,,, 为 多 年 时 间 序 列 中 i 月 地 表 温 度 的 
最 大 值 ;LS7,;, 为 多 年 时 间 序 列 中 i 月 地 表 温 度 的 最 
小 值 。 其 值 域 范 围 在 [0,1] 之 间 ,7CI 越 接近 于 0， 
该 月 越 干 旱 ,植被 生长 越 差 ,反之 亦 然 。 


TN 


TCI, = (2) 


2.1.3 PCI 指数 ”降水 指标 是 干旱 的 内 在 表现 ,长 
时 间 无 有 效 降水 或 降水 量 较 历年 同期 明显 偏 少时 通 
常会 发 生 旱灾 。 降 水 的 多 少 代表 干 湿 的 程度 ,干旱 
年 份 的 降水 量 与 正常 年 份 的 降水 量 往往 差异 较 大 。 
选取 与 VCI 及 TCI 类 似 的 基于 降水 时 间 序 列 变化 的 
PCI?) 作为 模型 参量 , 应 用 热带 测 两 任务 卫星 
(TRMM ) 数 据 参与 计算 , 即 


pare TRMM, - TRMM,,,, (3) 
"  TRMM pax — TRMM,,,, 

式 中 :PCL; 为 i 月 的 降水 状态 指数 ;TRMM; 为 i 月 的 
降水 量 ;TRMM,, 为 多 年 时 间 序 列 中 i 月 降水 量 的 最 
大 值 ;TRMM,;, 为 多 年 时 间 序 列 中 i 月 降水 量 的 最 小 
值 。 在 降水 较 历年 同期 均值 明显 偏 少 时 ,PCI 接近 
或 等 于 0, 表 明 越 干旱 ,反之 亦 然 。 

2.1.4 LC LC 对 干旱 的 发 生 也 有 一 定 的 影响 , 通 
HAN IAI) LC 类 型 土壤 含水 量 会 存在 差异 。 与 荡 地 
相 比 ,植被 能 够 改善 土壤 结构 ,增强 土壤 的 持 水 能 
力 。 和 森林 和 林地 对 土壤 水 分 具有 较 高 的 保持 力 , 草 
地 和 农田 则 只 能 利用 表层 土壤 的 水 分 , 当 出 现 长 时 
间 的 降水 短缺 或 高 温 蒸 散 , 便 会 因 土 壤 表 层 含水 量 
下 降 、 供 水 短缺 发 生 干旱 。 本 文 依据 不 同 LC 类 型 
对 于 水 分 的 保持 能 力 和 对 干旱 的 响应 程度 ,参考 夏 
SCHEMA MCD12Q1 数据 的 等 级 赋值 法 将 土地 重 
分 类 为 范 地 农耕 地 等 6 类 ( 表 3)。 由 于 LC 类 型 为 
年 值 数据 , 则 每 月 的 LC 类 型 对 应 同年 的 重 分 类 结 
果 。 

2.1.5 DEM 地形 因素 对 土壤 水 分 有 着 重要 影 
啊 。 高 海拔 地 区 的 太阳 辐射 能量 交 换 和 营 散 发 过 
程 与 低 海 拔 地 区 差异 较 大 ,会 对 干旱 发 生 细 节 产 生 
影响 ”] 。 对 于 地 形 因素 而 言 ,研究 区 域内 海拔 较 高 
的 地 区 发 生 干 旱灾 害 的 几率 较 小 ,反之 海拔 相对 较 
低 的 河谷 地 区 则 是 旱灾 的 高 发 区 。 由 于 描述 地 形变 
化 的 数字 高 程 模型 (DEM) .LC 与 YC1、7CI 以 及 PCI 


表 3 LC 重 分 类 
Tab.3 LC reclassification 
MCD12Q1 IGBP 分 类 编码 重 分 类 结果 
16 He Bea Tie Hb 1 荒地 
12 农田 14 农田 自然 植被 的 交错 带 2 农耕 地 
8 木 本 稀 树 草原 .9 稀 树 草原 、10 草地 3 草地 
1 常 绿 针 叶 林 、2 常 绿 阔 叶 林 、3 落叶 针 叶 林 、 4 林地 


4 落叶 阔 叶 林 ,5 混 生 林 、6 郁 闭 灌 从 7 FENNEMA 
13 城市 和 建筑 区 5 城镇 建设 用 地 
0 水 体 \11 永久 湿地 、15 冰雪 6 水 体 


KT > 作 HHF 


AA IVa FADT 
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的 量 纲 不 同 ,对 计算 结果 影响 较 大 , 故 对 二 者 采用 极 
差 标准 化 处 理 , 即 


NX 一 区 


, ij min 


a (4) 


式 中 :i,j 为 经 纬度 ;xy 为 (zj) 处 像 元 标准 化 后 的 像 
元 值 ;wj 为 (i,7) 处 像 元 值 ;xve 为 同期 影像 像 元 最 大 
值 ;xm 为 同期 影像 像 元 最 小 值 。 
2.2 模型 构建 

以 2001—2015 年 逐 月 对 应 的 VCI, TCI, PCI 与 
同期 标准 化 后 的 LC 类 型 和 DEM 数据 为 指标 ,利用 
ArcGIS 进行 空间 主 成 分 分 析 (SPCA)'" ,提取 累计 
贡献 率 大 于 85% 的 前 n 个 主 成 分 ,对 主 成 分 进行 加 
权 求 和 计算 ,得 到 研究 区 干旱 综合 监测 指数 (DC- 
MI) ,权重 系数 为 各 主 成 分 对 应 的 贡献 率 (图 2) , 即 


(5) 


式 中 :PC, 为 第 i 个 主 成 分 ;a 为 第 i 个 主 成 分 对 应 
的 贡献 率 。 
2.3 ”模型 验证 及 等 级 划分 
2.3.1 模型 验证 ”为 了 验证 构建 的 干旱 综合 监测 
模型 的 实用 性 与 准确 性 ,采用 土壤 相对 湿度 与 逐 月 
标准 化 降水 蒸 散 指数 (SPED 作为 验证 指标 , 通 
过 与 模型 计算 结果 的 对 比分 析 , 实 现 对 模型 的 定量 
评价 ,检测 其 监测 效果 。 

土壤 相对 湿度 验证 ,通过 研究 区 内 拉萨 站 等 4 
个 农业 气象 站 点 2001 一 2013 年 (数据 截止 至 2014 
年 ) 数 据 , 在 剔除 缺失 值 后 ,利用 逐 月 平均 10 cm £ 
壤 相 对 湿度 对 干旱 综合 监测 模型 进行 验证 。 农 业 气 
象 站 点 对 应 的 DCMI 取 站 点 周边 3 x3 区 域内 的 像 
元 均值 ,以 此 来 消除 个 别 像 元 值 异 常 变化 带 来 的 误 
差 , 最 后 利用 SPSS 软件 进行 相关 性 的 计算 分 析 。 研 


DCN = SosPC. 
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图 2 模型 构建 流程 
Fig.2 Model building process 
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图 3 研究 区 气象 站 点 分 布 


Fig.3 Meteorological site distribution in the study area 


究 区 农业 气象 站 点 分 布 如 图 3 所 示 , 相 关 性 验证 结 
果 如 表 4 所 示 。 

SPEI 指数 验证 ,通过 研究 区 内 28 个 气象 站 点 
2001—2015 年 逐 月 SPEI 指数 对 模型 结果 进行 验 
证 。 通 过 提取 农业 气象 站 点 周边 3 x3 区 域内 的 
DCMI 像 元 均值 与 SPEI 指数 进行 对 比 ,应 用 SPSS 
软件 对 数据 进行 相关 性 的 计算 分 析 。1 ~ 12 月 相关 
性 验证 结果 如 表 5 所 示 。 

结果 表明 DCMI 与 土壤 相对 湿度 及 SPEI 存在 
一 定 的 相关 性 , 且 相 关 性 达到 极 显著 的 水 平 。 农 业 
气象 站 点 土壤 相对 湿度 与 DCMI 相关 性 均 通 过 0.01 
的 显著 性 检验 , 且 相 关系 数 均 处 于 中 等 相关 ;1 ~ 12 
HZH SPEI 与 DCMI 相关 性 均 通 过 0. 01 的 显著 性 
检验 ,上 且 多 数 月 份 相关 系数 为 中 等 相关 。 因 此 ,DC- 
MI 能够 较 好 地 反映 区 域 土壤 相对 湿度 与 SPET 的 变 
化 ,干旱 综合 监测 模型 具有 较 好 的 适用 性 。 

2.3.2 等级 划分 通过 参考 植被 健康 指数 (VHI)、 
SPEI 干旱 等 级 划分 方法 “~” 以 及 我 国 气象 干旱 等 


表 4 DCMI 与 10 cm 土壤 相对 湿度 的 相关 性 
Tab.4 Correlation between DCMI and relative humidity 


of 10 cm soil 


站 点 日 喀 则 拉萨 泽 当 林芝 
相关 性 0. 503 * 0. 439 * 0.537 ** 0. 423 ** 


注 :* * 代表 通过 0.01( 双 侧 ) 显 著 性 检验 


#5 DCMI 5 SPEI 的 相关 性 
Tab.5 Correlation between DCMI and SPEI 


月 份 1 2 3 4 5 6 
相关 性 0.333" 0.321" 0.515 0.550 0.563“ 0.479™ 
月 份 7 8 9 10 11 12 
相关 性 0.552 0.599" 0.511 0.503 0.363 0.202 


注 :* * 代表 通过 0.01( 双 侧 ) 显 著 性 检验 


| 
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表 6 干旱 等 级 划分 
Tab.6 Drought classification 
序号 干旱 等 级 VHI SPEI DCMI 对 应 状况 

1 ， 极端 干旱 VHI<10 SPEIS -2.0 DCMI <0. 2 土壤 出 现 水 分 长 时 间 严 重 不 足 , 地 表 植 物 干枯 ET IKE 
和 生态 环境 造成 严重 影响 

2 ”重度 干旱 10<VHI<20 -2.0<SPEI< -1.5 0.2<DCMI<0.3 土壤 出 现 水 分 持续 严重 不 足 , 出 现 较 厚 的 干 土 层 ,植物 萎 薄 , 叶 
片 干 枯 , 果 实 脱落 ,对 农作物 和 生态 环境 造成 较 严 重 影响 

3 HEFE 20<VHI<30 -1.5<SPEI< -1.0 0.3<DCMI<0.4 ”降水 持续 较 常年 偏 少 ,土壤 表面 干燥 ,土壤 出 现 水 分 不 足 ,地表 
植物 叶片 白天 有 鞭 芒 现象 ,对 农作物 和 生态 环境 造成 一 定 影响 

4 REFE 30<VHI<40 -1.0<SPET< -0.5 0.4<DCMI<0.5 降水 较 常 年 偏 少 ,地 表 空 气 干燥 ,土壤 出 现 水 分 轻 度 不 足 , 对 农 
作物 有 轻微 影响 

5 JER  VHI>40 SPEI > -0.5 DCMI=0.5 降水 正常 或 较 常 年 偏 多 ,地 表 湿润 ,无 干旱 现象 


级 划分 标准 ,对 DCMI 进行 干旱 等 级 划分 ,划分 
结果 如 表 6 所 示 。 


3 西藏 东南 部 历史 干旱 识别 与 时 空格 
局 分 析 


为 了 进一步 研究 拉萨 .日喀则 等 5 市 的 干旱 时 
空 变 化 特征 及 发 展 趋势 ,首先 利用 所 构建 的 干旱 综 
合 监 测 模型 ,对 2001 一 2015 年 逐 月 的 DCMI 进行 识 
别 计算 ,对 研究 区 1 km 空间 分 辩 率 干旱 空间 信息 栅 
格 数据 集 进 行 重建 ,并 对 DCMI 计算 结果 进行 干旱 
等 级 划分 。 通 过 统计 分 析 各 级 干旱 的 发 生 频 率 ,从 
时 间 和 空间 两 个 维度 , 对 拉萨 、 日 喀 则 等 5 市 
2001—2015 年 干旱 时 空 演化 特征 进行 分 析 ,揭示 气 
候 变 化 背景 下 区 域 干旱 的 发 生 规律 与 发 展 趋势 。 
3.1 干旱 频率 统计 分 析 

干旱 频率 是 指 区 域内 发 生 干 旱 的 频繁 程度 1。 


文中 基于 GIS 空间 分 析 与 历年 逐 月 的 干旱 监测 结 
果 , 对 研究 区 干旱 频率 进行 空间 分 析 。 其 原理 是 当 
某 像 元 达到 干旱 阔 值 后 则 认定 此 像 元 发 生 干 旱 过 
程 ,赋值 为 1 ,其 余 像 元 赋值 为 0 ,最 终 获得 历年 逐 月 
的 干旱 统计 数据 。 

3.1.1 全 年 干旱 频率 分 析 全 年 干旱 频率 的 统计 
结果 ( 表 7) 表 明 ,2001 一 2015 年 拉萨 .日喀则 等 5 市 
平均 干旱 频率 为 13.22% ,87. 59% 的 区 域 干旱 频率 
小 于 20% 。 干旱 频率 在 20% ~30% 间 的 区 域 面积 
达到 11. 98% , A 0. 44% 的 区 域 干 旱 频 率 大 于 
30% 。 由 图 4 可 以 看 出 ,日 喀 则 市 全 年 干旱 频率 最 
高 , 且 以 轻 度 干旱 为 主 , 约 1/3 的 区 域 干旱 频率 大 于 
20% ,其 中 拉 孜 、 谢 通 门 、 定 日 等 县 的 部 分 区 域 干旱 
频率 均 超过 30% , 部 分 地 区 干旱 频率 其 至 超过 
40% ;拉萨 市 干旱 频率 较 高 区 域 主要 位 于 中 部 和 西 
南部 ,东部 和 北部 地 区 干旱 频率 相对 较 小 ;昌都 市 与 
林芝 市 大 部 分 区 域 各 级 干旱 频率 均 小 于 20% ,其 


表 7 研究 区 2001 一 2015 年 干旱 频率 分 布 
Tab.7 Drought frequency distribution of cities in the study area from 2001 to 2015 


干 时 等 级 频率 区 间 拉萨 市 日 喀 则 市 林芝 市 昌都 市 山南 市 

A% 面积 / km? 比例/ % 面积 / kim? 比例 / % 面积 / km” 比例/ % 面积 / km” 比例 / % 面积 / km” 比例 / % 

干旱 总 频率 <10 6 503 22.01 11 638 6.49 84574 73.96 39018 36.01 40949 51.84 
10 ~20 21227 71.86 108312 60.39 29695 25.97 69001 63.68 36 288 45.94 

20 ~30 1 808 6.12 57181 31.88 81 0.07 332 0.31 1743 2.21 

30 ~40 1 0.00 2 199 1.23 1 0.00 1 0.00 5 0.01 

>40 0 0.00 26 0.01 0 0.00 0 0.00 0 0.00 

轻 度 干旱 频率 <10 15677 53.07 40 423 22.54 104 924 91.76 83 692 77.24 57964 73.39 
10 ~20 13 840 46.85 136257 75.97 9 421 8.24 24658 22.76 20965 26.54 

20 ~30 22 0.07 2676 1.49 6 0.01 2 0.00 56 0.07 

中 度 干旱 频率 <10 29537 99.99 176367 98.33 114351 100 108 352 100 78 978 99.99 
10 ~20 2 0.01 2 988 1.67 0 0.00 0 0.00 7 0.01 

20 ~30 0 0.00 1 0.00 0 0.00 0 0.00 0 0.00 
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(c) 中 旱 频率 (qd) 重 旱 及 以 上 频率 
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图 4 研究 区 2001—2015 年 干旱 频率 分 布 
Fig.4 Drought frequency distribution in the study area from 2001 to 2015 
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5 研究 区 1~12 月 不 同 干旱 频率 面积 占 比 变化 


Fig.5 Changes in the proportion of different drought frequency in the study area from January to December 
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图 6 1~12 H DCMI 变化 趋势 空间 分 布 
Fig.6 Spatial distribution of DCMI trends from January to December 
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中 ,林芝 市 大 部 分 区 域 干旱 频率 小 于 10% ;山南 市 
中 偏 北部 区 域 干旱 发 生 频 率 明显 高 于 南部 区 域 , 约 
1/2 区 域 干旱 频率 小 于 10% 。 

从 不 同等 级 干旱 频率 来 看 ,研究 区 内 轻 旱 频率 

LTE 30% 以 下 ,其 中 99. 46% 的 地 区 轻 旱 频率 小 于 
20% ,日 喀 则 市 大 部 分 区 域 为 轻 旱 频 发 区 ;中 旱 频 率 
则 在 25% 以 下 ,其 中 99. 41% 的 区 域 中 旱 频率 小 于 
10% , 易 发 区 域 集 中 在 日 喀 则 市 的 部 分 地 区 ;2001 一 
2015 年 全 年 重度 及 以 上 干旱 频率 相对 较 小 , 均 在 
10% 以 下 ,主要 集中 于 日 喀 则 市 和 昌都 市 的 中 部 及 
东部 地 区 。 
3.1.2 月 际 和 干旱 频率 分 析 从 图 5 可 以 看 出 , 拉 
萨 .日喀则 等 5 市 月 际 干旱 频率 差异 较 大 。 整 体 
而 言 ,全 年 干旱 易 发 生 于 1 月 ,8 月 以 及 11 月 等 月 
份 ,局 部 地 区 干旱 易 发 月 份 存在 差异 。 其 中 ,1 月 、 
11 月 干旱 频率 在 20% ~ 30% 之 间 区 域 面 积分 别 
占 全 区 的 8.19% 与 6.99% 。 夏 季 7 月 .9 月 干旱 
频率 在 10% 以 下 的 区 域 面积 相对 较 大 ,表明 7 月 
和 9 月 不 易 发 生 干 旱灾 害 。8 月 虽然 为 雨季 ,但 由 
于 时 常 出 现 高 温 少雨 天 气 导 致 干旱 灾害 的 发 生 ， 
这 也 是 易 发 生 干 旱 区 域 面积 相对 于 7 月 .9 月 明显 
增 大 的 原因 。 

从 不 同 区 域 来 看 ,日 喀 则 市 旱情 易 发 生 于 每 年 
的 4 月 .6 月 、10 月 与 11 月 ,在 这 4 个 月 中 ,全 市 
42% 以 上 区 域 干旱 频率 超过 10% ;拉萨 市 旱情 则 易 
发 生 于 每 年 的 2 月 ,3 月 ,干旱 频率 10% 以 上 的 区 域 
达到 全 市 面积 的 45% 以 上 ;山南 市 除 9 月 外 ,其 余 
月 份 干 旱 频 率 均 为 均衡 ,全 市 20% ~39% 的 区 域 干 
旱 频 率 处 于 10% 以 上 ;林芝 市 旱情 易 发 月 份 则 主要 
集中 在 1 月 8 月 、11 月 与 12 月 ;昌都 市 旱情 易 发 月 
份 则 主要 集中 在 1 月 5 月 与 8 月 ,在 这 3 个 月 中 ， 
全 市 40% 以 上 区 域 干旱 频率 超过 10% 。 

3.2 ”旱情 变化 趋势 分 析 

为 了 定量 分 析 年 内 各 月 旱情 的 变化 趋势 ,利用 
一 元 线性 回归 分 析 方法 ,对 逐个 像 元 的 变化 趋势 进 
行 模拟 ,以 像 元 单元 的 变化 特征 来 综合 反映 区 域 时 
空格 局 变化 。 对 时 间 自 变量 及 DCMI 因 变 量 数据 ， 
采用 二 乘法 计算 所 有 像 元 的 DCMT 与 时 间 的 回归 和 斜 
率 , 即 


ni x DCMI, - = 


i=l i=l i=l 


Slope = 


式 中 :Slope 为 像 元 DCMI 回归 方程 的 斜率 ,i 为 月 序 
号 ,n 为 月 跨度 。 当 Slope >0 时 ,表示 DCMI 随时 间 
变化 呈 增 加 趋势 ; 当 Slope <0 时 ,表示 DCMI 随时 间 
变化 呈 下 降 趋 势 , 目 值 越 小 ,旱情 增加 趋势 越 明 显 。 

结果 表明 ,拉萨 .日 喀 则 等 5 市 年 内 旱情 变化 趋 
势 差 异性 较 大 (图 6)。 整 体 而 言 ,10 月 ~ 次 年 9 
月 ,旱情 加 剧 区 域 呈现 出 随 月 份 变 化 由 耕作 区 东部 
向 西部 逐渐 转移 的 趋势 。 其 中 ,10 月 ~ 次 年 3 月 ， 
山南 .林芝 .昌都 3 市 旱情 加 剧 趋势 较为 明显 ;4 ~5 
月 旱情 加 剧 区 域 面 积 达到 65% 以 上 ,加 剧 区 域 主 要 
集中 在 日 喀 则 .拉萨 .山南 的 部 分 地 区 ;6 ~9 月 , 拉 
萨 .日喀则 等 5 市 则 均 出 现 不 同 程度 的 旱情 加 剧 趋 
势 ,8 月 旱情 减弱 趋势 区 域 占 比 则 要 明显 高 于 其 他 
月 份 。 


4 结论 


本 文 以 地 处 西藏 自治 区 东南 部 的 拉萨 日 喀 
则 等 5 市 为 研究 区 , 以 2001 一 2015 年 MODIS, 
TRMM 遥感 数据 和 SRTM-DEM 为 数据 源 , 利用 
VCI TCI 和 PCI 等 模型 参量 ,采用 空间 主 成 分 分 析 
方法 构建 了 人 研究 区 干旱 综合 监测 模型 ,对 模型 的 
精度 和 可 靠 性 进行 检 核 验证 ,在 此 基础 上 ,对 
2001—2015 年 研究 区 逐 月 干旱 进行 识别 和 重建 ， 
分 析 了 2001 一 2015 年 研究 区 干旱 时 空格 局 及 变 
化 趋势 ,得 到 结论 如 下 : 

(1) 研究 区 域内 , DCMI 能 够 较 好 地 反映 区 域 
土壤 相对 湿度 与 SPET 的 变化 ,干旱 综合 监测 模型 具 
有 和 较 好 的 适用 性 。 

(2) 研究 区 全 年 干旱 频率 在 空间 分 布 上 呈现 
出 西高 东 低 的 特征 。 研 究 区 平均 干旱 发 生 频 率 为 
13.22% ,87.59% 的 区 域 干旱 频率 小 于 20%。 F 
旱 频率 在 20 ~ 30% 间 的 区 域 面积 占 11.97% , 约 
0.44% 的 区 域 干旱 频率 大 于 30% 。 从 不 同等 级 干 
旱 发 生 频 率 来 看 ,日 喀 则 市 为 轻 旱 、 中 旱 易 发 区 ， 
重 旱 易 发 区 则 集中 于 日 喀 则 市 和 昌都 市 的 中 部 及 
东部 地 区 。 

(3) 人 研究 区 月 际 干旱 频率 空间 格局 差异 较 大 。 
整体 而 言 ,全 年 干旱 易 发 生 于 1 月 .8 月 以 及 11 月 
等 月 份 ,局 部 地 区 干旱 易 发 月 份 存在 差异 。 

(4) 人 研究 区 年 内 旱情 变化 趋势 差异 性 较 大 。 整 
体 而 言 ,10 月 ~ 次 年 9 月 ,旱情 加 剧 趋势 呈现 出 随 
月 份 变化 由 耕作 区 东部 向 西部 逐渐 转移 的 趋势 。 
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Monitoring and analysis of historical drought in southeast Tibet based 


on multi-source data 
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Engineering Institute , Chengdu 610072 , Sichuan ,China; 4 Agriculture Research Institute , Tibet Academy of Agriculture and 
Animal Husbandry Sciences , Lhasa 850000 , Tibet , China ) 


Abstract: As one of the global natural disasters which take place frequently , drought has caused problems in soci- 
ety ,economy , and ecological environment. Based on the MODIS,TRMM and SRTM-DEM data from 2001 to 2015 as 
the data sources ,this paper takes the main farming area of Tibet Province ,China as the research area and construc- 
ted a regional drought comprehensive monitoring model by using the spatial principal component analysis method 
and adopting the vegetation state index ( VCI) , temperature condition index ( TCI) and precipitation state index 
(PCI). The accuracy and reliability of the model is verified. The model was then used to identify the monthly 
drought in the study area from 2001 to 2015, and the geospatial-temporal analysis method was used to study the 
characteristics and trends of drought changes in the study area. The results show that the drought comprehensive mo- 
nitoring index (DCMI) can better reflect the changes of regional soil relative humidity and standardized precipitati- 
on evapotranspiration index ( SPEI) ,and the comprehensive drought monitoring model has good applicability. The 
spatial distribution of annual drought frequency in the west of the study area was higher than that in the east, and the 
drought frequency in most areas was less than 20% ,and about 12.41% of the regions had a drought frequency more 
than 20% . From the perspective of the frequency of droughts in different grades, Shigatse City was a light drought 
and moderate drought area, while the severe drought areas were concentrated in the central and eastern parts of Shi- 
gatse City and Changdu City. The spatial pattern of monthly drought frequency in the study area was quite different. 
The drought in the whole year was prone to occur in January , August and November, and there were differences in 
the drought-prone months in some areas. The change trend of drought in the study area was quite different during a 
year. From October to September, the drought-intensified area showed a trend of gradually shifting from the east to 
the west of the cultivated area along the month. 
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